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The specific heat due to line nodes in the superconducting gap of YBa2Cu30 



7 



has been blurred up to now by magnetic terms of extrinsic origin, even for 
high quality crystals. We report the specific heat of a new single crystal 



grown in a non-corrosive BaZr03 crucible, for which paramagnetic terms are 
reduced to ~ 0.006% spin-1/2 per Cu atom. The contribution of line nodes 
shows up directly in the difference C(B,T) — C(0, T) at fixed temperatures 



(T < 5 K) as a function of the magnetic field parallel to the c-axis (B < 14 
T). These data illustrate the smooth crossover from C ex T 2 at low fields 
to C oc TB 1 ' 2 at high fields, and provide new values for gap parameters 



l I which are quantitatively consistent with tunneling spectroscopy and thermal 

conductivity in the framework of d x 2_ y 2 pairing symmetry. Data for B along 



the nodal and antinodal directions in the ab-plane are also provided. The 



in-plane anisotropy predicted in the clean limit is not observed. 
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I. INTRODUCTION 

Many experiments tend to establish that the symmetry of the order parameter in high 
temperature superconductors (HTS) is d x 2_ y 2, with a possible minor s-wave admixture |I],§. 
Whether they probe the amplitude or the phase of the order parameter, these measurements 
have been generally restricted to the surface of the samples. Alternatively, specific heat 
(C) experiments have been used to search for bulk evidence of a non-conventional gap, 
taking advantage of their sensitivity to the low energy excitations of a system. In a d x 2_ y 2 
superconductor, the gap A(/c) vanishes and changes sign on lines of nodes in fc-space. The 
finite slope of the gap at the nodes causes a linear increase of the density of states (DOS) at 
low energy, N(E) oc \E\. As the electronic part of C/T is proportional to the DOS averaged 
over an interval ~ ksT about the Fermi level, it follows that C e i ectron /T = aT at T -C T c in 
zero magnetic field. This contrasts with the exponential law C e i ectron /T oc T~ 2 ' 5 exp(—A /T) 
that would apply at very low temperature to a fully-gapped s-wave superconductor. 

In a magnetic field, the energy of carriers circulating around a vortex is shifted by 
Doppler effect. This allows them to be excited above the local gap on the Fermi surface 
near the nodes. Detailed calculations show that a contribution C vortex /T = A c B l l 2 arises 
in the T = limit j3]. In a conventional isotropic s-wave superconductor, no significant 
contribution would be expected at low temperature from such a mechanism, since a field on 
the order of B C 2 is needed to shift the energy by an amount comparable to the gap Ao- 

Two other mechanisms may contribute to the low temperature specific heat of super- 
conductors in a magnetic field: the Zeeman shift f||| and localized levels in vortex cores. 
H The Zeeman and Doppler contributions to C/T scale as B / B c2 and (B / B^) 1 ^ 2 , respec- 
tively, so that the latter dominates at low fields. Localized levels in vortex cores also add a 
contribution oc (B /B c2 ). For a low-T c superconductor with a large coherence length £, the 
mini-gaps ~ Ao//c_f£ between these levels are so small that the local DOS can be considered 
as continuous; the average DOS is proportional to the area occupied by vortex cores and 
entails C/T oc (B/B c2 ). This is the main field-induced contribution at low temperature 



for classic superconductors. For HTS, the localization in smaller vortex cores increases the 
separation between levels. Tunneling spectroscopy 0] tends to show that the mini-gaps are 
then on the order of ~ 0.3&bT c , which implies that the contribution of core levels should 
be negligible in the temperature range T < T c /20 investigated here. Furthermore, the very 
existence of localized core levels for a pure <i-wave superconductor is questionable, as they 
can leak along the nodal directions and mix with other states. 

Summarizing, the presence of line nodes in HTS is reflected by terms C(T <C T c , B = 
0)/T = aT and C(T = 0, B <C B c2 )/T = A C B 1 ^ 2 , whereas for a fully gapped, low tempera- 
ture superconductor C(T < T c , B = 0)/T oc exp(-A /T) « and C(T = 0, B < B c2 )/T oc 
B. However these simple criteria are not easily applied in practice because of the presence of 
additional contributions to the measured specific heat. In YE^CuaOy^, the main contri- 
butions arise from lattice vibrations, weakly interacting paramagnetic centers, and a linear 
term Cu n = 7(0)T of uncertain origin. All three are sample- dependent. In order to separate 



the (i-wave contribution of interest, the Stanford |8]|| and Berkeley |T0|JTl[] groups fitted mod- 
els to their data, and concluded that a <i-wave contribution was present. Numerical results 
were: a = 0.10 to 0.11 mJ/K 3 mol, A c = 0.88 to 0.91 mJ/K^T 1 / 2 mol according to Moler et 
al., |||S| who used both twinned and untwinned crystals, and a = 0.064 ± 0.02 mJ/K 3 mol, 
A c = 0.91 mJ/K 2 T 1//2 mol according to Wright et al., [|ll|] who used ceramic samples. In 
this paper, we report low temperature specific heat measurements for a twinned and fully 
oxidized YBa2Cu307 crystal in higher fields. The concentration of magnetic impurities for 
this new crystal is so low that it is possible to see for the first time in YE^CusOy-^ the con- 
tribution due to line nodes in the raw data prior to any correction or fit. Most uncertainties 
due to the background are avoided by looking at the difference between the specific heat in 
a magnetic field and that in zero field, C^n {B, T) = C(B, T) — C(0, T), which is electronic in 
origin. We obtain a = 0.21 mJ/K 3 mol±20% and A c = 1.34 mJ/K 2 T 1/2 mol±3%. Together 
with recent results from thermal conductivity, |L2| photoemission, |13| and tunneling spec- 



troscopy, |T4[ these new values allow a numerically consistent picture of the anisotropy of 



the order parameter for YBa2Cus07 to be constructed. A second issue that can be studied 



for the first time owing to the low magnetic contamination is the anisotropy of the vortex 
specific heat, which attracted recently theoretical interest. [[15] 



II. EXPERIMENTAL RESULTS 

1. Sample 

The sample used for the present measurement is a 18 mg YBa2Cu30"7.oo single crystal 
(code AE429G) flux-grown in a BaZrC>3 crucible. The latter non-reactive material allows 



a slow growth with little contamination. [16| Final chemical purity is better than 99.995%. 



T7fl Large crystals have to be screened because of possible CuO/BaCu0 2 flux inclusions. 



The effect of this second phase is potentially huge since only 0.9% BaCu02 in weight would 



double the heat capacity of the sample at 3 K, with an intricate field dependence. JT8[ The 
present crystal, about 2 x 1.5 x 1 mm 3 in size, is twinned, with its smallest dimension along 
the c-axis. Detwinning was not attempted, as it was found that standard procedures used for 
crystals of that size lead to a degradation of their properties in terms of transition width and 
concentration of magnetic centers. The sample was fully oxygenated in 100 bar O2, 330C, 
for 200 hours. The calorimetric transition midpoint, T c = 87.8 K, is typical of overdoped 
YBa 2 Cu30 7 .oo- The unextrapolated peak-to-peak amplitude of the C/T anomaly at T c is 
61 mJ/K 2 mol, i.e. 4.5% of the essentially phononic background (Fig. [[]). This amplitude is 
reproducible for crystals grown in BaZrC>3 (cf. 18 mg crystal AE195G |19| and 40 mg crystal 
AE276G |2(|). Fully oxidized samples are preferred for the present purpose for two reasons: 



clusters of oxygen vacancies effectively act as residual impurities, |H] introducing scattering 

that could modify drastically the shape of the DOS at E = in a <i-wave superconductor, 

22]1 and a correlation has been observed between the concentration of oxygen vacancies and 



the amplitude of the Schottky anomaly, [f2^] which tends to mask the d-wave contribution. 

2. Specific heat: method and tests. 

The low temperature specific heat was measured for 1.2 < T < 5K and < B < 14T 
using a thermal relaxation technique. [|23| For calibration purpose, we first measured in B = 



two silver samples of 6N purity with masses 9 and 21 mg. This allowed us to determine 
independently the heat capacity of the addenda (sapphire sample holder, contributing part 
of the phosphor bronze suspension wires, epoxy, carbon film heater/thermometer, silicone 
grease for sample mounting) and the specific heat of silver, which could be then compared 
with reference data. The heat capacity of the smaller Ag sample was on the same order of 
magnitude as that of the YBa2Cu3C>7 crystal in the middle of the temperature range. We 
verified that the measured specific heat of Ag after subtraction of the addenda did not change 
significantly with the field. Silver, which has a small nuclear magnetic moment, presents over 
copper standards the advantage of a negligible nuclear heat capacity in the temperature and 
field range investigated here. Figure |2] shows the measured low temperature specific heat of 
silver, together with reference data. [24| The insert of Fig. shows the Sommerfeld constant 



7 and Debye temperature G>d obtained in fields B = 0, 4, 8 and 14 T. The results are 
consistent with accepted values Bd = 226 K and 7 = 0.65 mJ/K 2 mol. |24| Note that for 
YBa 2 Cu 3 07, the Debye specific heat C oc O5 3 is 6.7 times smaller (per gram-atom), so that 
the determination of 7(-B) is more accurate. 

Specific heat data can be measured both with increasing and with decreasing tempera- 
ture. The former are obtained from the response to switching on the heater, the latter to 
switching off. The heat capacity is given in both cases by: 



[R(T(t)) - i?(T(oo))]/ 2 = C{T{t))^- + f m k(T')dT' (1) 

at jt 00) 



where R{T) is the resistance of the carbon film used both as a heater and a thermometer, 
I is the current through the film, set to a constant value at time t > 0, C(T) is the total 
heat capacity at temperature T, T(oo) is the final asymptotic temperature for a given 
value of / (measured typically after ten time constants), k(T) is the thermal conductivity 
of the supporting wires, measured separately when dT/dt = using various heater powers 
and various base temperatures. We use the approximations T(t) = (T n + T n+ i)/2 and 
dT/dt = (T n+1 — T n )/(t n+ i —i n ) with 50 ms time intervals between the n th and (n+ l) th data 
acquisition, therefore a single relaxation provides numerous independent C(T) points. These 



data tend to accumulate at the end of a relaxation until the error on (dT/dt) -1 diverges; we 
discard data with |T n+1 — T n \ < 4 mK. Such accumulations are seen in Fig. |2| and Fig. f|near 
T 2 = 2, 5, 11 and 22 K 2 . The up/down procedure allows one to check that measurements 
are independent of the heating rate (i.e., there is no 'V2 effect" in the sense of Ref. ||25| ). On 
a fine scale, there is a time dependence of the specific heat of YB^CuaOy at the highest 
fields and lowest temperatures. This is due to the hyperfine specific heat of Cu nuclei, which 
are weakly coupled to the electrons [ 27H . Data at T < 1.5 K in B = 14 T are discarded for 



this reason. 

3. Specific heat of the YBa2Cus07 crystal. 

Figure (| shows the low temperature specific heat of sample AE429G in different magnetic 
fields, shown separately for S||[001] along the c-axis (Fig. ^a), for S||[110] along the ab- 
diagonal where a gap node is expected (Fig. f|b), and for B along both the a- and 6-axes (the 
sample is twinned), i.e. along the antinodal direction (Fig. |]c). For the sake of simplicity, 
the latter direction is called _B||[100]. These are raw data, i.e. only the heat capacity of the 
sample holder is subtracted, otherwise no correction is performed for possible BaCuC^/CuO 
impurities originating from flux, [IT3] paramagnetic centers with S = 1/2 and/or S = 2, 



|2l| nuclear hyperfine contribution, |27| etc. In a first approximation the specific heat is 
essentially the sum of two dominant contributions. The C p h = j3T 3 lattice term shows up as 
a field-independent slope in the C/T versus T 2 plot, whereas the field-dependent electronic 
term C e \ = 7(-B)T shifts the curves parallel to each other. 

The Schottky effect of residual paramagnetic centers causes a low-temperature upturn 
in low fields, and a characteristic maximum of Csch/T at T max [K] = 0.415S[T] for S = 1/2 
spins. At high fields, the amplitude of the maximum decreases whereas its width increases, so 
that the Schottky contribution finally merges into the background. By fitting the expression 



Nz 2 e z 2gSfi B B eff 



Csch ~ (ez + i)2> z ~ & y 6 ' Be ff ~ V B ext + B int ( 2 ) 

we find iV = 0.12±0.02 mJ/K.gat, i.e. 0.005 to 0.007% spins-1/2 per Cu atom, and gS = 1.1 



to 1.15. |28| The amplitude iV and the value of gS are determined for B = 4 or 8 T along the 



afr-plane, in order to minimize the relative vortex contribution and to locate the maximum 
of the Schottky peak in a convenient temperature range. In this fit, which only serves to 
characterize the Schottky contribution, we have assumed that the other terms are (3T 3 for 
phonons and j(B)T for electrons (a more general analysis is given in Section 3). Imposing 
then the same value of N for B = or 0.5 T, we obtain an equivalent interaction field 
Ant ~ 0.2 to 0.5 T. The Schottky term is small relative to the vortex term for B\\c, but is 
significant for B^Lc at intermediate fields. In particular, it explains why the raw data at 14 
T do not lie clearly above those in 8 T in Fig. |]b and Fig. §c, in spite of the monotonically 
increasing vortex term. In any case, high fields improve the ratio of the vortex specific heat 
over the Schottky specific heat. 

Extrinsic magnetic contributions for sample AE429G are unusually small. This is not an 
artifact. Measurements performed in our laboratory on other samples, including ceramics 
p3| and detwinned single crystals from various sources have shown a more usual typical 
concentration of « 0.1% spin-1/2 centers per copper atom. Such "impurities" make it 
difficult to separate all contributions, in particular at low fields where magnetic interactions 
between spins cannot be neglected. 

The present measurements are compared with earlier work on YE^C^C^^ single 
crystal || in FigfJ. The same fields and the same presentation as for Fig. 3 of Ref |§ are 
used, i.e. a (3T 3 term is subtracted, he present data differ in the following respects: the 
measurement range is shifted by a factor of roughly two toward low temperature where 
the phonon contribution is smaller; the /3 coefficient is ~ 20% smaller (possibly owing to 
hardenng of Cu-0 chain at full oxygenation, but also to a smaller high-order electronic 
contribution); the magnetic Schottky term causing deviations from a horizontal line in non- 
zero field is suppressed by a factor of 15; the residual linear term 7(0) is ~ 30% smaller 
(but similar to that of sample Ul in Ref. ||). These features, which can be summarized by 
saying that all non-ci-wave terms contribute less, facilitate the analysisand give more weight 
to the d-wave specific heat for the present crystal. 



III. DISCUSSION 

1. Field along the c-axis. 

In order to get a premilinary insight into the final results, we plot again the data of 
Fig|4|(a) (5||c) in Fig.|](d), after having subtracted the 8 Tesla curve taken as a reference. 
The remaining specific heat qualitatively shows the evolution of the field-dependent elec- 
tronic contribution C e /T. Anticipating the discussion, we expect the high field curves to be 
free of the anomalous C oc T 2 term, hence the choice of the 8 Tesla curve as a reference. 
Higher fields would only introduce more scatter. This plot only involves smoothing of the 
reference curve (residuals are shown by the 8 T data set), but neither parameter fit nor 
correction of any kind. Fig.^(d) evidences the parallel shift of the curves in high fields, 
C c /T = f(B), and the progressive appearance of a positive slope at very low fields, ending 
into a C e /T oc T term in B = 0. The small upturn at low temperature in the low field 
curves, together with the faint maximum barely observable at intermediate fields, can be 
both explained by the presence of a small residual Schottky contribution. At the highest 
fields and lowest temperatures, data obtained from separate relaxations at different veloci- 
ties dT/dt do not join smoothly because of the weak coupling of electrons with Cu nuclei. 

The analysis of a difference C(T, B) — C(T, -B re f) circumvents the problem of modeling the 
field-independent background (lattice specific heat, residual "linear term", etc,). However, 
in the present status of knowledge a field -independent sub-linear or logarithmic term is no 
more justified tha a linear one, so that we avoid any fit of the background. The following 
discussion will only depend on differences. However, for a comparison with theory, B Te f = 
is more convenient as a reference than B Ie { = 8 T. Therefore we shall concentrate on 
C <m (T,B) = C(T,B)-C(T,0). 

Figure |^ shows Cda(B) at fixed temperatures T = 2, 3, and 4 K, prior to correction 
for the magnetic Schottky effect. It is immediately apparent that the vortex specific heat 
increases approximately with the square root of the field. This behavior can be readily 
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explained by existing theory, up to the deviations from the square root law. We briefly 
recall the essential results. 



Volovik et al. P, |29|J30|| first pointed out that in the mixed state of a superconductor with 
line nodes, supercurrents around a vortex core cause a Doppler shift of the quasiparticle 
excitation spectrum. If the superfluid velocity is v s , the quasiparticle excitation spectrum 
E(k) is shifted by k ■ v s . This shift has important effects around nodes, where its value is 
comparable to the width of the superconducting gap. The density of states at the Fermi 
level is strongly affected: 

N(0)=nJ-^Jd 2 r5(E(k,r)-k-v s (r)) (3) 

The average superfluid velocity depends on the inverse intervortex distance, (v s ) oc 1/R(B), 
so that the integral is proportional to R(B) oc 1/B 1 ^ 2 . For B 3> B c i, the number of vortices 
n is proportional to B, so that N(0) oc B 1 / 2 and C c \/T oc B 1 ! 2 at T — >• 0. Thus one has two 
regimes, depending on whether the thermal energy k-^T is large or small compared to the 
typical Doppler energy, one that is quadratic in T [C c \ = aT 2 ) in zero field at T ^C T c , and 
one that is linear in T {C e \ =A C TB 1 ' 2 ) at zero temperature and B <C B c2 - More precisely, 
for a weak-coupling superconductor with d x 2_ y 2 symmetry, Kiibert et al. |JT] and Vekhter 
et al. IH5U32I obtained: 



Ceuooi] ,8 i /2 B 1/2 14(2tt) 1 / 2 k B T 2 B c2 /a 2 1/2 TB c2 

7n T V { B c2 /a?> + 15 l A M B ' '"' T C B^ ^ {) 

gdjogi] 27C(3) k B T 3ln2 A g TB X J 2 

7nT vr2 A + 2tt kBTBa/a^-'TcBV*^ [ > 

where £(3) = 1.202... and a is a constant of order unity depending only on the vortex lattice 
geometry. The latter constant is defined differently in Ref. JHJ] and p!5|j32| ; we have chosen 



here the notation of Refs. |L2] and |3T| (see note |33|). These equations may be rewritten 
in terms of the Fermi velocity v-p to avoid any reference to the ill-defined upper critical field 
of HTS. In this form, they are no longer restricted to the weak-coupling case. Keeping only 
the leading terms, we have then: 



geUOOl] 4fl hv * Rl/2 , $ /2fc B T 

"^^ = ¥p^ + -'^T^ <<1 (6) 

CeUooi] 2703) fee d# 2 fc B T 

7n T 7T 2 A ' ^ F B 1 / 2 K ' 

These asymptotic formulas extrapolate to the same result at a crossover temperature 
T CTOBS (B) given by: 



:2tt) 3 / 2 An , B 



T cross (B) 27C{3)kB ( Bc2/a 2) ( 8 ) 

The ratio between the experimental parameters A c and a depends only on the Fermi velocity: 

A _ hrriT^oiC/TB 1 / 2 ) T CTOSS {B) = Anh 

a ~ UmB^{C/T*) BW 27C(3)$ / V F l ' 

The full function across the crossover regime has not been calculated analytically, but has 



the following scaling property: |30|j34 | 



^«^ 1/2 = F C)[001] (z) (10) 



7n T v B ' 
where x = T/T cross (B). [^3| The empirical interpolating function 



F cJO oi](o:) = B 1/2 a(l + a; 2 ) 1/2 (H) 

allows one to map the full specific heat function (Fig. |6|). Far above the crossover tempera- 
ture, the Doppler shift can be neglected, so that only the bulk term C/T oc T arising from 
the V-shape of the DOS survives. Conversely, far below T cross (_B), only the field- dependent 
plateau in the DOS at E — » is probed by thermal excitations, so that the Doppler term 
C/T oc B 1 / 2 dominates. In any case one has assumed T <^i T c and B <C B c2 , which is 
satisfied in the present data limited to T < 5 K and B < 14 T. As we shall see below, 
T cross / 'B 1 ' 2 = 6.4 K/T when B\\c, so that measurements of Cdit(T,B) from 2 to 4 K 
in fields from 0.16 to 14 T probe the region 0.085 < x < 1.6, which spans both the high 
field/low temperature limit and the crossover regime x ~ 1. This implies that both constants 
A c and a can be determined independently from the data. 

10 



Since we are interested in the difference C<jif, the relevant scaling function is F^{(x) 
defined by: 

C d jf, [ooi] _ , B 1/2 27C,(?>)k B _ B 1/2 

7nJ -Dc2 7T A -D c2 

and the corresponding interpolation function becomes: 

F* f>m] {x) = {^) l ' 2 a[(l + x 2 ) x ' 2 -x\ (13) 

Figure [7| shows the scaling plot Cda/TB 1 ^ 2 versus T/B 1 ^ 2 corresponding to Eq.(|i"2|). The 
data for B\\c, measured at fixed temperatures T = 2, 3 and 4 K, and now corrected for the 
Schottky contribution(in the main frame, not in the insert), collapse onto a single curve, 
thus supporting the existence of line nodes. This test also shows that the additional energy 
scale introduced by possible impurity scattering is negligible compared to the other scales 
determined by thermal smearing and Doppler shift. Strong impurity scattering would lead 
to a breakdown of the scaling property, with N(0) decreasing ultimately as |.Blog-B| rather 
than B x l 2 at low fields. [[H|] A comparison with numerical calculations (Fig. 4 of Ref. JTT] ) 
shows that our data lie essentially in the clean limit. This also implies that the observed 
residual linear term 7(0)T, with 7(0) ~ 15% of the normal-state Sommerfeld constant j n , is 
mostly not caused by impurity scattering. 

Figure § shows how parameters A c and a can be extracted directly from the data. We 
plot Cdif/T versus B 1 ! 2 at T = 2, 3 and 4K. At high fields, data fall on parallel lines. This 
is the region 1 C 1 where 

Ca/.poi] ws L596g 2^! _ 3 .288^T ee KB 1 ' 2 - aT (14) 



T lx ^~ i ~"~ Bl' 2 A 



c2 







The slope determines the leading field-dependent term with A c = 1.34 mJ/K 2 T 1/,2 mol ± 
3%, whereas the parallel shift with changing temperature in the high field limit determines 
the bulk term with a = 0.21 mJ/K 3 mol±20%. Their ratio yields avF = 1.0 x 10 7 cm/s. To- 
gether with the maximum gap width Aq = 20 meV determined by scanning tunneling spec- 



troscopy, [14j we further obtain 7 n = 15 mJ/K 2 mol and B c2 /a 2 = 310 T. These results agree 

11 



with the re-evaluation by Kiibert and Hirschfeld of earlier data characterized by various scat- 
tering rates. |H] The crossover temperature is given by T cross / ' B 1 ! 2 = A c /a = 6.4 K/T 1 / 2 . 
The parameters A c and a are substantially larger than those determined earlier and re- 
called in the Introduction. This may be due in part to the use of a crystal that is overdoped 
rather than optimally doped. Scattering by impurities is probably lower. ]3T| In addition, 



the methods used by different authors have different sensitivities to various sources of ex- 
perimental errors. The present results allow one to construct a scenario in the framework of 
d x 2_ y 2 pairing symmetry that is consistent with other types of experiments such as tunneling 
spectroscopy and thermal conductivity. Following Chiao et al. [[Tj| we assume a cylindrical 
Fermi surface; the A c parameter can be expressed in terms of the slope v 2 = dA/dp± of 
the gap on the Fermi surface in the direction perpendicular to the line of nodes: 

8fc B Vmol O /-, r\ 

A c = TR— I (15) 

3/i$o /2 d V2 
This follows from Eq. |, with 7n = (ir 2 /3)klN(0)V moh JV(0) = m*/(nh 2 d) where d = 
c/2 = 0.584 nm is the average interlayer distance, V^oi — 104.6 cm 3 /mol the molar volume, 
A(0) = A o cos(20), and |dA(0)/d0|^, =7r /4 = hk F v 2 = 2A at the node. Eq. (15) yields 
v 2 /a = 1.4 x 10 6 cm/s. In the same approximation, a is given by 

18C(3)fc| VW 1 , 1R . 

a = 72 3 (16) 

nn d v 2 vp 

which yields v 2 vf = 1.4 x 10 13 (cm/s) 2 . Up to this point, the geometrical parameter a is 
undetermined. The recent analysis of thermal conductivity in the universal low temperature 
regime by Chiao et al. [12| determined v-p/v 2 = 14 ± 20% in YBa 2 Cu30~7. Combining this 



result with ours one obtains v F = 1.4 x 10 7 cm/s, v 2 = 1.0 x 10 6 cm/s, a = 0.70, B c2 = 150 
T, corresponding to a coherence length £ = 15 A. The geometrical parameter is of order one 
as anticipated, and the upper critical field is consistent with previous estimations. p5|j36[] 
The Fermi velocity is smaller than determinations based on angle-resolved photoemission 
spectroscopy (ARPES), v F = 2.5 x 10 7 cm/s. |TB| The difference appears to be real, even 



taking into account the large uncertainty margin on this parameter (Table I), and may be 
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due to a contribution from chain states with a low Fermi velocity in our overdoped sample. 
In this respect it is interesting to note that YBa 2 Cu40g, which contains completely filled 
double chains, presents a large negative curvature in the plot of C/T versus T 2 in B = 
at T < 5 K. |23| If this curvature is attributed to a bulk d-wave aT 2 term, then a is still 
three times larger (per chain) for YBa 2 Cu 4 08 than for YBa 2 Cu 3 07, and consequently the 
product t>2fF i s three times smaller, showing a potential source of variations of vp. The 
Fermi wave number is kp = 2Ao/frt>2 — 0.6 A = 3/4 7r/oi (ai = 3.85 A is the lattice 
constant). The latter result is consistent with the quasi-2D hole band centered on the S- 



point (k = [n/ai, n/b, 0]) found by band structure calculations |37|] and by ARPES, [13 



which show that plane states cross the Fermi surface near kp = 0.69 to 0.75 A along 
the nodal direction and kp = 0.58 to 0.64 A along the antinodal directions. The overall 
consistency provides strong support to the quasi-2D, d^.^-wave model. In particular using 
a normalized gap function A(0)/Ao similar to that of Bi2Sr2CaCu20s does not lead to any 
contradiction. 

Based on the previous data, one can further evaluate the Fermi temperature Tp = 
hk-pv-p i '2k-Q = (A /kB){vp/v2) — 3300 K and the mass enhancement factor m*/m e = 
hkp/m e vp = 5. Note that the latter ratio depends only on the product of experimental 
quantities Aq«. It is interesting to note that the Bose-Einstein condensation temperature 



of such heavy quasi-2D hole pairs, if preformed above T c , should occur at [38 



2irh 2 nd I 



Tbe = ~, „, 7r P2 ,,2 = 100 K (17) 

k B m* ln( 2fcBTB g md ) 



where T = 5.3 is the anisotropy ratio of the upper critical field. p9| , ^0| For the pair density, 
we used the value n = 1.0 x 10 21 cm 3 obtained by scaling the amplitude of the A anomaly of 
the specific heat of 4 He at 2.18 K, for which the boson density is n/, = 2.2 x 10 22 cm 3 , 
to the A-anomaly of YBa2Cu3C>7 at T c . The fact that Tbe is close to T c may be a 
coincidence for YBa2Cu30y, since the penetration depth calculated from these parame- 
ters exceeds experimental values by a factor of about two. Additionally, the ratio |T2 



Ao/Ep = 2/(7t/cfO — (^/^f) — 1/14 shows that only a small part of the carriers near the 
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Fermi surface are paired. 

Finally the initial assumption that vortex core levels do not contribute significantly to 

the specific heat at low temperature can be tested a posteriori. Self-consistent calculations 

of the electronic structure of a vortex line in NbSe2 have shown that the lowest energy 

excitation depends on temperature, with a value at T <C T c of the order of E1/2 — A 2 /Ep. 

4l| Applying the same criterion to YBa2Cu307 leads to Ei/ 2 /kB — 17 K, about half the 



estimation based on STS, @ but large enough to be neglected in the specific heat below 4 
K. 

The above determinations are summarized in Table |. In order to evaluate uncertainty 
margins, we included the error on experimental parameters A c , a, and v-p/v 2 , and neglected 
that on Ao- 

2. .B 1 1 aft-plane 

The main motivation of the measurements for I?_l_c was the prediction of a variation 
of j(T = 0, B, <p) with the azimutal angle <p between the field and the a-axis. [[HJ The 
qualitative physics underlying this effect may be understood by noticing that when a field 
is applied along the antinodal [100] direction, all four nodes contribute equally to the DOS, 
with an amplitude proportional to 4cos(7r/4), whereas for B along the nodal direction [110], 
the Doppler effect vanishes for those particles which travel parallel to the field, so that only 
two nodes contribute with an amplitude proportional to 2cos(0). Quantitative estimations 
depend on details of the theory, so that a 30% variation appears as an upper limit for a clean, 
tetragonal crystal with a 2D band. Such fourfold oscillations of j(T = 0, B, 0) versus 0, if 
observed, would provide a robust test of d-wave symmetry. Moler et al. || did not observe 
any such oscillation within experimental resolution. A second motivation is the experimental 
verification of the general features of the vortex specific heat predicted theoretically for I?_Lc 



HP1 



Data analysis for I?_Lc is more involved than for B\\c (compare Fig. |9] and Fig. |5]). Even 
for the present crystal, the residual Schottky contribution due to paramagnetic centers is 
no longer negligible compared to the reduced vortex contribution. We face for the 5J_c 

14 



configuration the same problem as the Stanford |8||| and Berkeley [|T0|JTTf] groups did for the 
B\\c configuration. As data for f?_l_c are expected to lie in the crossover regime (the crossover 
temperature scales with (B/T) 1 ! 2 ), a fit to a C =A a _\ ) TB 1 l 2 law is not justified beforehand. 
We resort then to a scaling plot which is more generally valid, first using raw data (Fig. 
|TU|a; note that the scaling variable is B 1 ! 2 jT in this plot, so that unlike Fig. [7] high fields 
are on the right). At this stage, data do not collapse on a common curve. In the same figure 
we plot the Schottky specific heat of spin-1/2 paramagnetic centers, with a concentration 
equivalent to 0.006% of the Cu atoms. The position of the peaks corresponds to the maxima 
of the raw data, confirming that the paramagnetic centers have a spin 5* = 1/2. Note that 
the Schottky correction is negligible in B = 14 T. In the scaling plot of Fig. |T0"lb, we show 
the data after having subtracted the Schottky contribution. The parameters of the Schottky 
anomaly were adjusted manually to obtain the best collapse onto a common curve. In this 
operation, the most efficient parameter is the concentration of paramagnetic centers; only 
minor improvements are obtained by refining the gyromagnetic ratio and the interaction 
field. Data define a wide high field region (T/B 1 / 2 < 2 K/T 1 / 2 ) where C = A^TB 1 / 2 , 
A &h = 0.18 mJ/K 2 T 1/2 mol ±10%, therefore A c /A ah = 7.4. The latter determination does not 
depend critically on the subtraction of the Schottky anomaly; it suffices to consider only 14 
T data, _B_l_c, where the Schottky contribution can be neglected, and compare with 14 T 
data, 511c (Fig. | and Fig. |). 

The ratio A c / ' A^ is larger than the square root of the anisotropy ratio at T c , T 1 ! 2 = 



2.30 ± 0.07. This is expected, at least qualitatively. We recall the results of Ref. 1 15,32 



for the specific heat with the field parallel to the planes, assuming moderate anisotropy, 

cylindrical Fermi surface and Ak = Aocos(20): 

C d ,2,,,,, BIT Nl/9 27C(3U B T lir . TiTB^) 1 / 2 , s 

^ = '-/'W + ^~t + ■■" B|l[1101 ' T$- « 1 <«) 



C el 2 , BIT sl/9 28^/2 k B T, 9/ B c2 /a 2 ^ /9 ... , T(TB c2 ) l/2 , . 



C„ _ 27((3) k B T , 31n2 A„ B/T „...„„. „,,,„„, T(TB. 



c2 



,1/2 



7n r- ^ a ■ 4, fe r B ^ + -- B|l|110| - B|l|1001 -^i^ >>1 (20) 
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As before, we have restored the convention for the geometrical parameter a of Ref. [31]; see 



Note [3"3"||. In the high temperature regime, one recovers the bulk isotropic C oc T 2 behavior, 



irrespective of the field direction. At T = 0, the expected ratio between the c-axis and 
nodal specific heat is C(0,fi||[001])/C(0,5||[110]) = 2T 1 / 2 = 4.6, that between the c-axis 
and antinodal specific heat is C(0, B\\ [001])/C(0, B\\ [100]) = (2/T) 1 ' 2 = 3.3. The theoretical 
parameter-free curves given by the second-order approximations, Eq. |T^-|20|, are plotted in 
Fig. |i~0b. The order of magnitude is correct, but Eq. |T^-|20| first fail to reproduce the wide 



plateau for B 1 ' 2 /T > 0.5 T 1 ' 2 /K, second tend to overestimate the vortex specific heat for 
large values of 5 1//2 /T, and third predict an in-plane anisotropy for large values of B 1 / 2 /T 
that is not observed. There is no way to improve the agreement by changing the spin and 
the amplitude of the Schottky correction, which is negligible for large values of B/T. 

One is tempted to conclude that the anisotropy ratio at low temperature is larger than 
the value V = 5.3 found at T c (an example of such behavior is given by 2H-NbSe2 for which 



(£ c2 ,ab/£ c2 , c ) = 3.4 at T ^ whereas (dB c2:Sih /dT)/(dB c2 JdT) = 2.7 at T c ), g2j but this 
does not remedy the situation. A larger anisotropy would indeed decrease the coefficient of 
the vortex term A a ^, and therefore would shift down the theoretical curves; however this 
would leave the coefficient of the bulk term a unchanged and, remembering that the zero 
field specific heat is subtracted in Cdif, this bulk term would cause a positive slope for all 
values of B 1 / 2 jT plotted in Fig. pT0|b. Therefore, the plateau would not be reproduced. 

The experimental points show that the difference between the in-field and the zero- 
field specific heat contains essentially a small vortex term A^TB 1 ! 2 for T jB 1 ! 2 < 2 K/T 1//2 , 
apparently without any — aT 2 contribution from the bulk term. Leaving aside the disturbing 
scenario a = 0, one concludes that the application of a field 5±c leaves the aT 2 term 
unchanged at low temperature, so that it cancels in Cdif. In terms of density of states, a 
small V-shape dip subsists at E F in the middle of a field-induced plateau that moves up with 
B 1 ' 2 . Within this picture, the sharp crossover seen experimentally near T/B 1 ' 2 = 2 K/T 1 ' 2 
corresponds to the edge of the dip, not to A^/a = 0.86 K/T 1 / 2 . Such a fine structure of 



the DOS could indicate that the anisotropic London theory leading to Eq. |18| - [20|, and which 
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is valid for moderate anisotropy, does not take into account all aspects of vortex physics for 
B±c. 

The absence of in-plane anisotropy is a second puzzling feature. It is already apparent 
in the raw data, Fig. [|b and Fig. |]c, and more precisely seen in the scaling plot of Fig. [L0]b. 



where the points for B along the antinodal direction hardly differ from those for B along 
the nodal direction. In order to perform a more quantitative analysis, we focus on the 
B — 14 T data for two reasons: first the Schottky correction is negligible, and second the 
vortex contribution follows a simple high field <i-wave law, as just discussed. By fitting 
the parameters of a model C/T = 'y(B) + j3T 2 to the data in B = 14 T (Table |I|), we find 
7(.B|| [100])— j(B\\ [110]) = 0.04±0.06 mJ/K 2 mol. The residuals show no structure suggesting 
any missing term. The sign is correct but the amplitude is an order of magnitude smaller than 
the d-wave prediction, A^ 1 ' 2 ^)- 1 ' 2 - (4F)" 1 / 2 ] = 0.45 mJ/K 2 mol at 14 T in the clean 
limit. The 8 T data confirm the 14 T result. In order to make this comparison as significant 
as possible, we took care not to add or remove any addendum mass (in particular the 
adhesive) when rotating the sample at room temperature. The experimental reproducibility 
is documented by the data in B = 0, which were measured separately in each one of the 
three sample positions [001], [100] and [110] along the magnet axis (Table [TJ). 

In order to reconcile the lack of any significant variation of 7(0) with the existence of 



d-wave pairing, various arguments may be invoked. [32| First, the 2D model of the DOS may 
be too crude. Numerical estimations have shown that the amplitude of the oscillations of 
j(T = 0, B, 0) with is reduced in the 3D case. |32| Second, the full 30% effect develops only 
at T = 0; distinct crossover effects decrease the in-plane anisotropy at any finite temperature, 
|15| but this reduction is taken into account in the model curves of Fig. pTOjb. Third, 
orthorhombicity shifts the sharp minimum of j(T = 0, B, = 7r/4) slightly off the nodal 
direction. |32[ Together with twinning, this replaces the single minimum by a double dip 



centered on = tt/A with reduced depth. Scattering may smooth the oscillations; |43| 
we believe that this is not the main cause here since scattering would also have led to a 
breakdown of the 7 ex B lj/2 dependence for B\\c. Finally, the implications of the in-plane 
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anisotropy of the penetration depth in YBa2Cu307, A a > At,, [@] were not addressed. In 



short, YBa 2 Cu 3 07 does not appear to be suitable for a quantitative test of d-wave pairing 
symmetry based on afe-plane anisotropy, owing to its orthorhombicity. 

IV. CONCLUSION 

The low contamination of the present crystal by paramagnetic impurities has allowed 
to observe directly a field-induced contribution to the specific heat of YE^CusOy that can 
be accounted for by the theoretical treatment of a o?-wave superconductor with a quasi- 
2D band at the Fermi level. The parameters inferred from this experiment have a typical 
±20% accuracy as they depend critically on the value of the bulk T 2 -term of the zero- 
field specific heat. Together with the results of tunneling spectroscopy, photoemission and 
thermal conductivity, these results allow an overdetermined and consistent set of microscopic 
parameters to be established. 

The bulk confirmation of ci-wave properties provided by specific heat is important, as 
stressed by previous studies on this subject. However experiment and theory still indi- 
cate some unsolved puzzles. First recall that a fully gapped s-wave component, if present, 
would not give rise to any measurable contribution at T <C A /£;b, so that the present 
experiments cannot provide any information on this topic. It has been pointed out that 
non- linear variations of j(B) with the field also occur in s-wave superconductors, both from 
an experimental |E|,(|^,f|(J and theoretical [47 , 48] point of view. The present experiments 



probably cannot distinguish between C oc B 1 ' 2 and C oc B 0A1 as proposed in Ref. |48 



The expected afe-plane anisotropy could not be observed. As the present results leave little 
hope of obtaining a decisive answer using YBa2Cu307 as a working material, tetragonal 
compounds should rather be investigated, preferably with a Ba-free composition in order 
to avoid the presence of BaCuC>2 impurities. The field dependence of the specific heat for 
Bl.c remains to be explained in detail. Finally, it has been pointed out that the relatively 
large local DOS that should be present between the vortex cores in the Doppler model has 



not been observed up to now by tunneling spectroscopy in YBa2Cu307, notwithstanding 



several attempts | 4"5|j . Therefore, in spite of considerable progress, a full understanding of 
the mixed state of YBa2Cus07 has yet to be reached. Bulk investigations of unconventional 
pairing in other HTS, if they can be prepared with the same degree of purity, would be most 
informative. 

V. APPENDIX 



In a previous article, referred to as BR, |50[ we have estimated by a different method 
the vortex specific heat of another crystal (AE37G), also grown in BaZrC>3, with a larger 
Schottky anomaly corresponding to 0.03% spin-1/2 per Cu atom. The field-independent 
terms were canceled by considering the anisotropic component C(B\\c) — C(i?_Lc) rather 
than the difference C(B) — C(0). It was assumed that anisotropy enters only through 
the ratio -B C 2,ab/-Bc2,c = 5.3, an estimation that now has to be revised. BR obtained 
A c = 1.8 mJ/K^T^mol, and, noticing that [C(B\\c) - C(BLc)]/TB 1 ' 2 was not constant, 
concluded that data were taken in the crossover regime. 

Now that the specific heat for l?_Lc has been measured separately, we can check 
these conclusions. The present finding is that the specific heat for B^Lc is much smaller 
than expected. Rather than C(B\\c)/C(B±c) oc (-B c2iab /S c2jC ) 1/2 = 2.3 as used in BR, 
C(B\\c)/C(B-Lc) = 7.4 over a wide field range. Therefore C(B\\c) — C(i?_Lc) does not con- 
tain 1 — r -1 / 2 = 57% of C vortex (-B||c), but 1 — 1/7.4 = 87%. Our previous result, based on a 
measured anisotropic component A c — A a j = 1.0 mJ/K 2 T 1 ' 2 mol at T — > 0, has to be rescaled 
to A c = 1.2 mJ/K 2 T 1 / 2 mol. The difference with the present result, A c = 1.34 mJ/K 2 T 1//2 mol, 
is essentially due to the way data were extrapolated to T = in Fig. 3 of BR. 

The second point is the crossover temperature that separates the low T, high B regime 
from the high T, low B regime. In a plot of [C(B\\c) - C^B^c^/TB 1 / 2 versus T/B 1 / 2 , the 
crossover is the value of T/B 1 / 2 at which [C(.B||c) — C(5_Lc)]/Ti? 1//2 extrapolates linearly to 
zero, starting from its initial value at T/B 1 / 2 = (see e.g. Fig. |7|). In view of the smallness of 
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the vortex contribution for £?_l_c, one has approximately C(B\\c)—C(BJ_c) ~ C(B\\c) — C(0). 
Fig. 2 of BR, redrawn on a linear scale, shows that the crossover estimated in this way is 
A c /a = 5 to 7 K/T 1 / 2 for B\\c, consistent with the result for the present crystal AE429G, 
A c /a = 6.4 K/T 1 / 2 . In the high field region < T / B l l 2 <C A c /a, thus below the crossover 
temperature, C(B\\c)/TB 1 ^ 2 is constant but [C(B||c) — C(0)]/Ti? 1/ ' 2 must decrease linearly 
with T/5 1//2 , a point that was not made clear in BR. 
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After completion of this work, we became aware of recent developments in theory, both 
in the Doppler (Volovik) approach |51] and in self-consistent calculations of the electronic 



structure of a d x 2_ y 2 vortex p2"|j . Although the shape of the envelope is preserved, the density 
of states at the Fermi level in the latter approach is always zero, on a very fine energy scale, 



irrespective of the field [|53|-|55|j. This is not in contradiction with the present results in the 



1-4 K range, which are broadened by thermal smearing. 
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FIGURES 
FIG. 1. Total specific heat of the present crystal (AE429G) and another crystal also grown in 

BaZrO"3 (AE195G) near T c . The critical temperature is defined here as the point with the largest 

negative slope of C/T versus T. 

FIG. 2. Specific heat C/T versus T 2 of the silver test sample. Full line: standard reference 



data. [24 1 Insert: fitted Debye temperature and Sommerfeld constant in different magnetic fields. 



FIG. 3. Comparison of the present data(series (b) on the left side 1.2 - 2.4K). As for Fig.|| of 
Ref. H|, the specific heat is plotted per mol units as (C — f3T 3 )/T vs T, (3 = 0.392mJ/molK for 
the data of Ref. ||, and (3 = 0.305mJ/molK for the present work. The magnetic field is applied 
perpendicular to the planes. The same subset of fields is shown. 

FIG. 4. Low temperature specific heat C/T of YBa2Cu307 (AE429G) for different fields, raw 
data, (a), £||[001], from bottom to top B = 0, 0.16, 0.5, 1, 2, 3, 4, 6, 8, 10, 12, 14 T. (b), B||[110], 
B = 0, 0.5, 1, 2, 4, 8, 14 T. (c), B||[100] and [010], B = 0, 8, 14 T. (d), same data as (a), after 
subtraction of the B = 8T curve 

FIG. 5. Specific heat C/T versus the field B (B\\c) at fixed temperatures, raw data, showing a 
nearly square-root law. 

FIG. 6. Model calculation of the low temperature specific heat versus temperature and field 
(-B||c) in the presence of line of nodes, showing the high- field, low temperature regions where 
C/T oc B 1 ' 2 and the high temperature, low field regions where C/T oc T. The parameters 
correspond to those for the present crystal AE429G. 

FIG. 7. d-wave scaling plot of the specific heat difference C^ii/TB 1 ' 2 versus T/B 1 ' 2 (B\\c). A 
Schottky contribution has been subtracted. Full line: interpolated scaling function. Insert: the 
same plot without correction for the Schottky contribution. 
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FIG. 8. Plot of the specific heat difference C^a/T versus B 1 ' 2 (B\\c) at fixed temperatures, 
allowing one to extract the parameters A c and a. A Schottky contribution has been subtracted. 
Full lines: interpolated scaling function. Insert: the same plot without correction for the Schottky 
contribution. 

FIG. 9. Plot of the specific heat difference C^/T versus B (S||[110]) at fixed temperatures, 
raw data. 

FIG. 10. (a), same raw data as for Fig. 9 (.B||[110]) in a scaling plot C^/TB 1 / 2 versus B l / 2 /T. 
The full lines show the estimated Schottky contribution at different temperatures, (b), remaining 
vortex contribution after having removed the Schottky anomaly, both for Sj|[110] and S||[100] 
Dotted and dashed lines: anisotropic d-wave model (see text). 
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TABLES 



sample 

Vmol = 13V^at 

M mo \ = 13M gat 



r gg 

A c 

A ah 

a: 

v 2 /a 

VpV 2 

avp 
a 2 vp/v2 



T IB 1 / 2 

± crossover / L> 



A 

In 



B c2 /a 2 
m* /m e 
v F /v 2 IT 



Vf 
V2 
B C 2 

T F 



AE429G 

104.6 

666.2 

87.8 ± 0.05 

5.3 ± 0.3 

1.34 ± 3 % 

0.18 ± 10 % 

0.21 ± 20 % 

1.42 x 10 6 ± 3 % 

1.4 x 10 13 ± 20 % 



1.0 x 10 7 ± 24 % 



6.9 ± 27 % 
6.4 ± 24 % 

20 

15 ± 20 % 

310 ± 11 % 

5.1 ± 20 % 

14 ± 20 % 
0.70 ± 23 % 

1.4 x 10 7 ± 53 % 

1.0 x 10 6 ± 27 % 

150 ± 68 % 

15 ± 30 % 
0.61 ± 27 % 
3300 ± 20 % 



cm 3 /mol 

g/mol 

K 

mJ/K^ 1 / 2 !!^ 

mJ/K^ 1 / 2 !!^ 

mJ/K 3 mol 

cm/s 

(cm/s) 2 

cm/s 

K/T l/2 

meV 

mJ/K 2 mol 

T 



cm/s 
cm/s 
T 

A 
A- 1 

K 



2N 



TABLE I. YBa2Cu3C>7 parameters obtained from the field-induced specific heat. Rows below 
Ao make use of the gap value from tunneling spectroscopy. @ Rows below v-p/v2 make use of 
this ratio from thermal conductivity. |1^| See text for definitions. The effect of experimental 
uncertainties on A c , a and v-p/v2 is indicated; the error on Aq is neglected. 
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-B(T) orientation j(B) [mJ/K 2 mol] 0d (K) r.m.s. error [mJ/K 2 mol] number of data 

357 
647 
305 
514 
301 
210 
723 

TABLE II. Results of fits of the high field specific heat for l?_Lc: coefficient of the linear term 
'y(B) and Debye temperature 0d- The Schottky correction is fixed and corresponds here to 0.007% 
spin-1/2 per Cu atom with g = 2.2 and B nit = 0.5 T. The data in B = with 3 orientations with 
respect to the magnet axis are shown only to document the reproducibility, since the bulk C oc T 2 
d-wme term is not included in the fit (attempts to include it yield values that are strongly correlated 
with i?i nt and 0d and do not improve significantly the fit). 
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